Four-body continuum-discretized coupled-channels calculations 
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The development of a continuum-bins scheme of discretization for three-body projectiles is nec- 
essary for studies of reactions of Borromean nuclei such as ^He within the continuum-discretized 
coupled-channels approach. Such a procedure, for constructing bin states on selected continuum 
energy intervals, is formulated and applied for the first time to reactions of a three-body projectile. 
The continuum representation uses the eigenchannel expansion of the three-body S-matrix. The 
method is applied to the challenging case of the ®He + ^"^^Pb reaction at 22 MeV, where an accurate 
treatment of both the Coulomb and the nuclear interactions with the target is necessary. 

PACS numbers: 21.45.-v,24.50.+g,25.60.-t,24.10.Eq,27.20.+n 



Rapid experimental developments have made studies 
of the scattering and reactions of very weakly-bound 
systems possible in many branches of physics. Among 
these are the Efimov states observed in ultra-cold ce- 
sium trimers [l| and molecular triatomic systems, such 
as LiHe2 halo nuclei Q, systems undergoing two- 
proton radioactivity Q, and also Bose-Einstein conden- 
sates and ultra-cold dilute gas studies. The need for novel 
and quantitative theoretical descriptions of the dynam- 
ics of such systems, which must include a priori a reahs- 
tic treatment of their continuum spectra, is common to 
many problems. Here we discuss continuum effects in the 
nuclear physics context of reactions of weakly-bound halo 
nuclei 0, a The study of these systems is of special 
importance in nuclear physics in order to learn proper- 
ties of the nucleon-nucleon interaction in the low density 
limit and investigate the role of pairing off the stability 
line. Many of these exotic nuclei play a significant role in 
astrophysical processes. Although the method is general, 
in this Rapid Communication, we will discuss reactions 
of the Borromean halo nucleus ^He, i.e. the bound three- 
body (^He+n+n) system with (i) no bound excited states 
and (ii) no bound binary sub-systems. 

Due to their weak binding, halo nuclei are readily 
excited to the continuum (broken-up) by the differen- 
tial forces exerted on the constituents through the nu- 
clear and Coulomb interactions with a target. Theo- 
retically, the explicit treatment of these breakup chan- 
nels is difficult since the states involved are infinite in 
number and are not square-integrable. Hence a robust 
discretization and truncation scheme, to provide a finite 
and normalizable basis to represent the continuum, is 
required. The continuum-discretized coupled-channels 
(CDCC) method, as formulated for two-body projec- 
tiles 0, uses bin states to represent the continuum. 
Here, the continuum spectrum is truncated at a maxi- 
mum excitation energy £max and the model space, from 



the breakup threshold to £max, is then divided into inter- 
vals, where the number and positions of the intervals can 
depend on the properties (e.g. resonant or non-resonant) 
of the continuum of the system. For each such interval, 
or energy bin, a representative square-integrable state 
is constructed as a linear superposition of the two-body 
scattering states in the interval. The method has been 
enormously successful in the description of elastic and 
breakup observables in reactions involving weakly-bound 
two-body projectiles [1, [13, EH and has been recently 
extended to include core excitation [l^ . 

Of interest here are reactions of three-body projec- 
tiles. To date, published discretized three-body pro- 
jectile calculations have used a pseudostate (PS) basis 
15, 16], accounting reasonably well for existing 
elastic scattering data. Based on our own experience 
[ibI of calculations using PS bases, convergence prob- 
lems were found for reactions where the Coulomb inter- 
action plays an important role. The alternative contin- 
uum treatment using the energy bin technique within the 
CDCC has yet to be formulated and the present work 
provides this theoretical development. 

To describe the three-body ground and excited contin- 
uum states of the projectile we make use of a hyperspher- 
ical harmonics (HH) expansion basis This involves 
use of the one radial and five angular hyperspherical co- 
ordinates, p, a, X, y, obtained from the normalized Jacobi 
coordinates x, y of the three bodies [H, [l3|. The quan- 
tum number set, (3, that defines each three-body chan- 
nel [l3| are the hypermomentum K, the orbital angular 
momenta Ix and ly in coordinates x and y, their total 
I = Ix + ly, the total spin Sx of the particles associated 
with coordinate x, and the intermediate summed angular 
momentum jg^jj — 1 + Sx- If the spin of the third particle, 
/, is assumed fixed then the total angular momentum is 
j = Jab + ^ ^ith projections fi. Note that, for each con- 
tinuum energy e and total angular momentum j, there 
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will be as many independent solutions of the three-body 
scattering problem, as the number of outgoing channels 
(3 considered. These solutions can be chosen as the in- 
coming channels /3', but any orthogonal combination of 
these could be equally valid. 

Based on these total angular momentum eigenstates, 
3^/3 j>(^^) where Vt = (a, x, y), the bin wave functions 
are defined as 

^'^%ix,y) = J2R^:;,)ip)yp,,{n), (i) 

where the label n includes reference to the energy inter- 
val of the bin [ki , K2] , as well as to the set of quantum 
numbers /?'. The functions R^^jip) in Eq. ^ are the 
associated hyperradial wave functions, 

Rnfjip) = R\n'^,K2]!3'l3jiP) (2) 

where k = y/2m\e\/h is the momentum associated to 
the continuum energy e, Rf3f3'j{Kp) are the continuum 
hyperradial wave functions with their scattering 

phase shift, and fpij{n) is a weight function with Npij its 
normalization constant. Note that for each j and three- 
body energy bin we must construct wave functions for all 
allowed incoming channels /?'. Further, as is explicit in 
Eq. Ill) , for each n we must also construct R^pj [p) for all 
allowed outgoing channels (i. 

It follows that to include a large number of (3' chan- 
nels is a severe computational challenge, and that it is 
desirable to establish a hierarchy of the continuum states 
according to their importance to the reaction dynamics. 
In so doing, we may be able to describe scattering ob- 
servables using only a selected set of states, the number 
of them depending on the reaction under study. To this 
end, we make use of the eigenstates of the multi-channel 
three-body S-matrix fl8|, or eigenchannels (EC), as fol- 
lows, (i) For each j and continuum energy e, the S-matrix 
in the (3 basis is diagonalized to obtain its EC, enumer- 
ated by 7, their corresponding eigenvalues exp[2i(5.yj (k)] 
and eigenphases 5^j{K). (ii) The magnitudes of these 
eigenphases are used to order the EC. We will show that 
those EC with largest phase shifts are the most strongly 
coupled in the reaction dynamics, and thus a hierarchy of 
states can be established by such an ordering. This leads 
to the possibility of a truncation in the number riec of 
EC included and testing of the convergence with respect 
to this number. 

Here we apply this methodology to the scattering of 
^He, treated as a three-body system of an inert a parti- 
cle core and two valence neutrons. A notable property of 
^He is that none of its binary sub-systems bind, while the 
three-body system has a single bound state with bind- 
ing energy of 0.973 MeV and total angular momentum 
j'^ = 0+. Its low-lying continuum spectrum is dominated 
by a narrow j'^ = 2+ resonance, 0.825 MeV above thresh- 
old. We describe the three-body a+n+n system with the 



same structure model as used in Ref. [l5|. The three- 
body Hamiltonian includes two-body interactions plus an 
effective three-body potential. For a given value of the 
maximum hypermomentum used, -fCmax, the parameters 
of the latter are adjusted to reproduce the ground state 
separation energy and matter radius (for j = 0+) and 
the resonance energy (for j = 2+) j'15]. Calculations were 
performed using the codes FACE [19] and STURMXX [20| . 
The maximum hypermomentum in the CDCC reaction 
calculation was -ftTmax = 8, which provides converged re- 
sults for the elastic scattering of ^He from a heavy target 
[l5| . The number of channels (3 was 15 (0+), 26 (1~) and 
46 (2+). The calculated ^He ground state has binding 
energy of 0.953 MeV and a single-particle density with 
root mean squared (rms) radius 2.46 fm, assuming an a 
particle rms radius of 1.47 fm. 

The EC basis states have properties that are useful to 
describe collisions. The upper panels of Figs. [1] and [2] 
show the Coulomb B(E1) and B{E2) transition proba- 
bilities, respectively, as a function of the ^He excitation 
energy above threshold. The curves show the total B(E1) 
and B{E2) strengths (thick solid) and the contributions 
from the first four EC. Since the ^He(g.s.) has = 0+, 
B{E1) measures the electric dipole strength to 1^ states, 
and B{E2) the quadrupole strength to 2+ states. For 
each EC, the 1~ and 2+ eigenphases are also shown in 
the lower panels of the figures. It has been shown (see for 
instance Ref. j2l|) that the position of the peak of the 
B{E1) distribution depends on the maximum hypermo- 
mentum chosen. In particular, for -fsTmax = 10 the max- 
imum appears around 2 MeV, and decreases to about 1 
MeV, for i^max = 20, for which convergence of this ob- 
servable is achieved. For i^max — 8, which is the value 
used in this work, the maximum is around 2 MeV, and 
hence this observable is not converged. Nevertheless, as 
we have shown in Ref. the scattering calculations 
are converged with this value, provided that the param- 
eters of the three-body potential are adjusted such that 
the ^He system has the same binding energy and rms 
matter radius for the chosen i^max- In principle, the re- 
action calculations could be done with a higher iiTmax 
value, although this would turn these calculations com- 
putationally very demanding. 

It is evident from Figs. [T] and [2] that there is energy 
localization for the EC contributions to the EX strength. 
Also evident is that those EC with eigenphases of the 
largest magnitude make the maximum contribution to 
the B{EX) at low excitation energy. This is a very ap- 
pealing feature for scattering studies where, at moderate 
collision energies, only continuum states up to a few MeV 
play a significant role in the reaction process. So, for ex- 
ample, if excitation energies up to « 8 MeV are strongly 
populated in a particular reaction, the first three or four 
EC should suffice to describe the relevant part of the con- 
tinuum. Of course, in nuclear collisions near and above 
the Coulomb barrier, nuclear forces (not included in Figs. 
[U and [21) will play an important role. We show below 
that, with nuclear interactions present, only a few EC 
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Figure 1: (Color online) B{E1) distribution (upper panel) 
and eigenphase shifts (lower panel) for the first few = 1~ 
EC of ''He. The thick solid line in the upper panel is the total 
B{E1) distribution. 
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Figure 2: (Color online) B(E2) distribution (upper panel) 
and eigenphase shifts (lower panel) for the first few j'^ = 2^ 
EC of ®He. The thick solid line in the upper panel is the total 
B{E2) distribution. 



are needed for converged scattering observables. 

It is also interesting to identify the dominant quantum 
numbers of the EC with the largest contributions to these 
observables. For B[E1) the most important component 
of EC 1 and 2 has Inn = Ix = 0, Sx = 0, or a dineutron 
configuration. This dominance, could explain that previ- 
ous two-body {a+2n) model calculations of ^He reactions 
work to a limited extent, although it is well known that 
such simple models overestimate continuum coupling ef- 
fects without parameter adjustments [22, 23). For B{E2) 
the most relevant structure is the 2+ resonance that cor- 
responds to the EC 1. In this case we have a com- 
bination of three different configurations: a dineutron 
{Ix = Q,Sx = 0), a cigar-like structure {Ix = 2,Sx = 0), 
and a spin-triplet structure {Ix = 1, = 1). 



We now compute the states of Eq. ([T]) and use these as 
the basis states for a four-body CDCC calculation using 
the techniques of Ref. There, a multipole expansion 
of the coupling potentials was developed for a three-body 
projectile plus target system. The weight function f-fj{K) 
(with 7 instead of /?') used in the superposition of states 
when constructing the bins Eq. ^ was taken to be unity 
for the non-resonant (0+, 1~, and 2+) continuum and 
sm{5^j{K)) for the resonant (2+) continuum. The latter 
provides an improved description of the resonant charac- 
ter. 

Four-body CDCC calculations are carried out for the 
^He + ^"^^Pb reaction at 22 MeV, for which experimen- 
tal data are available from Ref. [2J| . This is one of the 
more demanding examples, due to (a) the importance of 
long-range Coulomb couplings arising from the heavy tar- 
get and (b) the incident energy being near the Coulomb 
barrier and Coulomb and nuclear forces both playing a 
significant role. 

Implicit in the CDCC is that the actual continuum 
can be truncated at a maximum excitation energy Emax, 
whose choice is dependent on the specific reaction and 
projectile energy. The convergence of the model space 
calculation must then be studied with respect to £max 
and the number of bins assumed (nbin)- For our four- 
body case we have, in addition, to study convergence with 
respect to riec, the number of EC included. The present 
calculations include ^He states with j'^ — 0+,!^ and 
2+ and projectile-target interaction multipole couplings 
with order Q = 0,1,2. Coulomb and nuclear potentials 
are included. The nuclear interactions with the target 
used parametrized optical potentials. For n-|-^°*Pb and 
Q,_^208p|3 ^j^g potentials were from Refs. j^Bl and [2i|, 
respectively. 

The coupled equations describing the projectile-target 
motion were solved with the code FRESCO [27||, the cou- 
pling form factors being read from external files. Partial 
waves up to J = 150 were included and the solutions were 
matched to their asymptotic forms at radius Rm — 200 
fm. In order to check the convergence with respect to 
the basis size, for the different observables shown in this 
work, we have used five different sets that are summa- 
rized in Table HI Sets I, II, and III have the same max- 
imum energy £max and number of EC rioc, so they are 
used to test the convergence with respect to the number 
of bins considered for each EC and j'^. Set IV is like set 
II but including the fifth EC, so it will be used to study 
the convergence with n^c. Finally, set V is like set II but 
including an additional bin with e=8-9 MeV for all EC 
and j"^ , providing information on the convergence with 

respect to Emax- 

First, we analyze the elastic scattering. Fig. [3] com- 
pares the calculated and experimental (circles) elastic 
differential cross section angular distributions (ratio to 
Rutherford). The dashed line is the one channel cal- 
culation, in which continuum states are omitted. The 
solid (set I), the dotted (set II), dot-dashed (set IV) and 
double dot-dashed (set V) lines are the full CDCC cal- 
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Set 


£max (MeV) 


Woe 


nbin(0+) 


nbm(l~) 


nbin(2+) 


N 


I 


8 


4 


6 


9 


6 


85 


11 


8 


4 


9 


12 


9 


121 


III 


8 


4 


12 


15 


12 


157 


IV 


8 


5 (1-4,5) 


(9,6) 


(12,9) 


(9,6) 


142 


V 


9 (0-8,8-9) 


4 


(9,1) 


(12,1) 


(9,1) 


133 



Table I: Different sets of states used in this work for the four- 
body CDCC calculations. N is the total number of states 
considered (including the ground state). The parentheses are 
used to specify the two different bin schemes adopted depend- 
ing on the EC (Set IV) or on the energy (Set V). 
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Figure 3: (Color online) Elastic differential cross section (ratio 
to Rutherford) in the center of mass frame for the ''He+^^^^Pb 
reaction at 22 MeV. See text for details. The experimental 
data are taken from Ref. [24] . 

culations with different choices of the model space. It is 
remarkable that all these calculations are almost indistin- 
guishable, showing that the elastic angular distribution 
is converged with set I. Also, this full four-body CDCC 
calculation, in contrast to the one channel calculation, 
describes the data fairly well. Set III is not shown since 
the convergence with the number of bins is already seen 
with sets I and II. 

Solution of the CDCC coupled equations for the 
projectile-target S-matrix also provides predictions for 
breakup observables. First, Fig. |4] shows the breakup 
differential cross section angular distribution, summed 
over the excitation energy. The solid (set I), dotted (set 
II), dashed (set III), dot-dashed (set IV) and double dot- 
dashed (set V) lines are the full CDCC calculations with 
different choices of the model space, as in Fig [31 For 
this breakup observable the rate of convergence is slightly 
slower than for the elastic cross section, although the con- 
vergence reached with set II is fairly good. 

Second, Fig. [5] shows the breakup cross section dis- 
tribution, as a function of ^He excitation energy above 
breakup threshold, for each j'^ and for each EC. This 
distribution is calculated by dividing the breakup cross 
section to each bin by its energy width. We use here set 
II that gives a converged result for the breakup angular 
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Figure 4: (Color online) Breakup differential cross section 
angular distribution in the center of mass frame for the 
''He+^°^Pb reaction at 22 MeV. See text for details. 



distribution. As was anticipated, in all j'^, the contri- 
bution of each EC decreases with decreasing magnitude 
of its eigenphase. In particular, the contribution of the 
fourth EC is already small, justifying the use of a trun- 
cated calculation and Uec ~ 4. We have also checked that 
the contribution of the fifth EC (comparing sets II and 
IV) changes the total breakup cross section by less than 
2%. This figure also shows the importance of the dipole 
and resonant states. These give the largest contribution 
to the breakup cross section. 

Collisions of loosely bound three-body projectiles with 
a target can be studied in the continuum-discretized 
coupled-channels framework. The continuum-bin dis- 
cretization procedure has been formulated and used for 
the first time for a three-body projectile, requiring the 
superposition of three-body scattering states. For each 
three-body total angular momentum j and energy e, scat- 
tering states are calculated for ingoing boundary condi- 
tions with quantum numbers (3. The multi-channel S- 
matrix in this /? basis, determined from the asymptotics 
of these states, is diagonalized to determine the eigen- 
channels (EC) that are ordered according to the magni- 
tude of their associated eigenphases. We have shown that 
those EC with the largest eigenphases are most strongly 
coupled in the collision, suggesting a hierarchy of contin- 
uum states and an associated truncation scheme. This 
EC hierarchy allows practical CDCC calculations that 
extend the energy binning procedure to four-body reac- 
tions induced by three-body projectiles. The formalism 
has been applied to the '^He+ 2°®Pb reaction at 22 MeV. 
The convergence of the results for elastic and breakup 
cross sections have been presented. Good agreement with 
existing elastic data for this reaction has been obtained. 
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Figure 5: (Color online) Breakup cross section distribution for the ®He+^'^*Pb reaction at 22 MeV with respect to ®He excitation 
energy (over the breakup threshold) for each j'^ and EC included in the calculation. Basis set II was used. 
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